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Swift observations of SAX J1808.4-3658: monitoring the return to 
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ABSTRACT 



The transient accreting millisecond pulsar SAX J1808.4-3658 has shown several 
outbursts to date but the transition from outburst to quiescence has never been in- 
vestigated in detail. Thanks to the Swift observing flexibility, we monitored for the 

^ '. first time the decay to quiescence during the 2005 outburst. At variance with other 

transients, wide luminosity variations are observed. In addition, close to quiescence, 
SAX J1808.4-3658 seems to switch between two different states. We interpret them 

Trt ', in terms of the accretion states accessible to a magnetized, fast rotating neutron star. 



Subject headings: accretion, accretion disks — star: individual (SAX J1808.4-3658) 
r" ! — stars: neutron 

'NT 

f^ ■ 1. Introduction 

O 

oo . 

O ■ SAX Jl 808.4-3658 (SAX J1808 in the following) was the first neutron star transient during 

whose outbursts coherent pulsations were detected (Wijnands & van der Klis 1998). This source, 

K> . together with six other low mass transients (see Wijnands 2005 for a review), form a separate (sub- 

H ! )class within neutron star transients (see Campana et al. 1998a for a review). Distinctive properties, 

besides coherent pulsations, are weaker outburst peak luminosities (~ (1 — 5) x 10^^ erg s"^), very 

low mass companions (mass functions < 2 x 10^^ M©), short orbital periods (Porb ~ 4.3 hr), very 

faint quiescent luminosities ( ^ 5 x 10^^ erg s"^) and absence of a soft X-ray spectral component 

in quiescence. 
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A further difference with respect to 'classical' neutron star transients concerns the return to 
quiescence following an outburst. The best light curve to date of 'classical' neutron star transients 
is represented by BeppoSAX observations of Aql X-1 (Campana et al. 1998b). Below a lumi- 
nosity level of ~ 10'^^ erg s~^ Aql X-1 turned off with an exponential decay with an e— folding 
time of ~ Id. The source then remained quiescent for at least one month following the 1998 
outburst (Campana et al. 1998b) and for five months following the 2000 outburst (Rutledge et al. 
2001; Campana & Stella 2003). The return to quiescence for transient pulsating neutron stars has 
been monitored for SAX J1808 (Wijnands et al. 2003; Wijnands 2005) and for XTE J1751-305 
(Markwardt et al. 2002) and Swift J1756.9-2508 (Krimm et al. 2007). The behaviour of the latter 
two sources was somewhat similar to Aql X- 1 with a fast decay, even though in the case of XTE 
J1751-305 there was still activity 15 d after the transition to quiescence had occurred (according to 
RXTE/PCA, corresponding to a 0.5-10 keV unabsorbed luminosity of ~ 10^^ erg s~^ for a source 
distance of 7 kpc). In case of SAX Jl 808 the decay was completely different from that observed in 
Aql X-1. The most striking feature is the strong erratic variability by a factor of ^ 30. This erratic 
behaviour was not concentrated in the first few weeks after turn off but lasted for months (e.g. 
Wjinands et al. 2004). The most likely explanation for this behaviour comes from the similarities 
with dwarf novae where a similar behaviour has been observed (Kato et al. 2004; Patterson et al. 
2002) and to a lesser extent in black hole transients (Kuulkers, Howell & van Paradijs 1996). In 
the case of dwarf novae this has been interpreted as due to a combination of short orbital period 
and low mass ratio (Osaki & Meyer 2002; Truss et al. 2002; Campana et al. 2008, in preparation). 

Here we report on a monitoring campaign carried out with Swift, following the 2005 outburst. 
At variance with previous monitoring with RXTE/PCA we are able to follow the source down to 
very low levels, a factor of ~ 50 fainter than before, close to quiescence. In Sect. 2 we discuss the 
data and their analysis. Sect. 3 contains our discussion and conclusions. 



2. Data analysis 

Swift carried out 23 observations of SAX J1808 between Jun 17, 2005 and Oct 28, 2005. 
XRT collected a total of 4899 s data in Window Timing (WT) mode and 72922 s data in Photon 
Counting mode (see Table 1). In WT mode a ID image is obtained reading data compressing along 
the central 200 pixels in a single raw. PC mode produces standard 2D images (for more details see 
Hill et al. 2004). WT data were mainly collected during the early stages of the outburst when the 
source was brighter. PC data collected during this period are piled-up. 

The XRT data were processed with standard procedures (xrtpipeline ver. 0.11.6 within FTOOLS 
in the Heasoft package ver. 6.4), filtering, screening, and grade selection criteria (Burrows et al. 
2005). For the WT data, we extracted source events in a square region with a side of 20 pix- 



Table 1: Observation log. 



Obs. ID. Date WT exp. time PC exp. time Count rate Counts 





(2006) 


(s) 


00030034001 


17Jun 


119 


00030034002 


20Jun 





00030034003 


23Jun 


761 


00030075001 


29Jun 


1043 


00030034005 


07 Jul 


866 


00030034006 


13 Jul 


733 


00030034010 


02 Aug 





00030034011 


05 Aug 





00030034012 


13 Aug 


1028 


00030034013 


21 Aug 





00030034014 


28 Aug 


68 


00030034015 


30 Aug 





00030034016 


31 Aug 





00030034017 


02 Sep 





00030034018 


14 Sep 





00030034019 


25 Sep 


123 


00030075020 


30 Sep 


42 


00030034021 


04 Oct 





00030034020 


11 Oct 





00030034022 


12 Oct 





00030034023 


16 Oct 





00030034024 


20 Oct 





00030034025 


28 Oct 


116 



(s) 
787 
1061 
206 




394 
416 
1041 
2486 
1675 
1928 
2550 
2838 
2084 
9011 
7329 
4922 
7757 
2986 
4689 
6742 
4736 
7284 



(c s-i) 

4.46 

1.49 

3.01 

4.09 

8.14 

3.10 
8.88E-03 
6.31E-03 

1.84 
9.05E-03 

1.36 
6.26E-02 
9.87E-03 
5.48E-03 
1.50E-03 
6.47E-03 
4.89E-03 
1.74E-03 
4.89E-03 
1.64E-03 
1.83E-03 
<3.49E-03 
<3.00E-03 



3510 (P) 

1581 (P) 

620 (P) 

4396 

7157 

1221 (P) 

4 [A] 

7 [A] 

4574 (P) 

15 [A] 

2622 (P) 

160 
28 [A] 

11 [A] 

14 [B] 
47 [A] 
21 [A] 

12 [B] 

15 [A] 
7[B] 
12 [B] 

16 
22 



Labels in the last column indicate: (P) piled-up photon counting data; [A] state A observation; [B] 
state B observation. 



els. Ancillary response files generated with xrtmkarf and accounted for different extraction 
regions, vignetting, and point-spread function (PSF) corrections. In PC mode we extracted data 
with variable extraction regions depending on source strength, ranging from annular region with 
inner (outer) radius of 10 (40) pixels when the source was very bright, down to 10 pixel circular 
region when the source was very faint in order to increase the source signal vs. the background. 
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Fig. 1.— SAX J 1808 light curve observed with Swift XRT during the 2005 outburst. 



2.1. Outburst light curve 



The first observations caught SAX J 1808 during the latest stages of the outburst phase at a 
count rate level (corrected for pile-up) of more than 10 counts s^^. Starting from MJD 53564 the 
source decreased its count rate by a factor of ~ 600 in 20 days (unfortunately no observations 
were taken during this period due to the occurrence of number of gamma-ray bursts, see Fig. 1). 
Two observations within 3 days showed the source at a level of 0.01 counts s^^ well above the 
quiescent level. After this sharp decay, SAX J 1808 entered a "flaring" behaviour (Wijnands et al. 
2003; Wijnands 2005). Our data evidence that this activity lasted for about one month; after that 
the source returned to its quiescent state with some comparatively small scale variability (factor of 
~ 5) superimposed (see also below). 



-5 



10 
1 



^S^^^ 



i Outburst 




0.5 1 2 

Energy (keV) 



Fig. 2. — XRT spectra of SAX Jl 808.4-3658. The topmost spectrum is a representative of the 
outburst state. The middle spectrum refers to state A and the lowest one to state B (see text for 
more details). All spectra were fitted with an absorbed black body plus power law model. 



2.2. Spectral analysis: bright end 

We fit the entire data set witli an absorbed black body plus power law model, keeping the same 
value of the column density for all the spectra. The best column density is (1.3±0.1) x 10^^ cm^^ 
(90% confidence level, modelled with TBABS). In Table 2 we report the results of our spectral 
modelling. From this table and Figs. 2 & 3 it is clear that during the outburst decay and flaring 
period the power law photon index (F ~ 2.3) and the black body temperature (/cT ~ 0.65 keV) 
remained almost constant. This behaviour occurs in the (unabsorbed) 0.3-10 keV luminosity range 
of 7 X lO^'' — 10^^ erg s~^ (i.e. a factor of ~ 15), assuming a source distance of 3.5 kpc (Galloway 
& Cummings 2006). 



2.3. Spectral analysis: faint end 

At lower luminosities the number of collected photons is very small and we have to stack 
together different observations. One important consideration comes from the observation that 8 
out of 13 observations found the source with a count rate in the narrow interval 5 — 10 c ks^^, 
whereas the other 5 all lie within 0.8 — 2c ks^^. The mean of the first group is 6.3 ± 1.1 c ks^^ 
and of the second is 1.3 ± 0.4 c ks~^. This is rather peculiar and hints for the existence of two 
different states in the deep faint end tail of the outburst. Due to the small number of photons these 
spectra can be easily fit with single power law models. For state A we derive a photon index of 
r = 2.7 ± 0.3 and a mean unabsorbed 0.3-10 keV luminosity of 5.0 x 10^^ erg s~^, for state B we 
have r = 1.7 ± 0.6 and a luminosity of 1.5 x 10^^ erg s~^. These values have been derived binning 
the data to 5 photons per energy bin and applying the Churazov weighting in the fitting procedure 
(see Fig. 2). This indicates that state B is not yet the true quiescent state, characterized by a single 
power law spectrum with index r = 1.93]']q 29 and a source luminosity of (5. 2±0.1) x 10^^ erg s^^ 
(Heinke et al. 2007), i.e. a factor of ~ 3 smaller that state B luminosity. Assuming this power law 
component as a stable component (at least in photon index) present in state A and B spectra, one 
can investigate if a black body component is present. Fitting state A spectrum with a fixed power 
law photon index (with F = 1.93) and a free black body component, we derive kT = 0.20 ± 0.05 
keV and a radius of R = 1.3^q 5 km. The improvement over the simple power law fit is at a level 
of 2.2 (T by means of an F-test. For state B we are not able to constrain the black body component. 
However, if we assume the same black body temperature as in state A (i.e. 0.2 keV) we can derive 
an upper limit on its radius of i? < 0.3 km, indicating that if present and at the same temperature 
it must be smaller in size. 



Table 2: Spectral fits. 

Date BB Temperature BB Radius PL index Luminosity'^ 

(keV) (km) (erg s^^) 

2005-06-17^ 0.72+°:{| 2.3+^:|^ 2.24to:n 1-0 x 10^6 

2005-06-20^ O.SQtoll l-S+ag 2.51+°;^^ 2.4 x lO^^ 

2005-06-23^ 0.59+^;}^ 2.4+^;^ 2.26+[J:J^ 7.2 x lO^^ 

2005-06-29^ 0.78+^;}^ L0+°:^ 2.37+[J:}^ 3.3 x lO^s 

2005-07-07^ 0.44+°:| 2.2+2-2 2.41+°;°^ 5.7 x lO^s 

2005-07-13^ 0.57l[J;J^ 2.2+[J| 2.18+[J:J^ 4.2 x lO^s 

2005-08-13^ 0.68+°:?^ 0.9+°| 2.22+°;^^ 2.0 x lO^^ 

2005-08-28^^ O.rStofr O-Stpi ^-Q^-om 1-4 x 10^5 

State A^ 0.20+°:°^ L3+°j 1.93 fixed 5.0 x 10^2 

State B'i 0.20 fixed < 0.3 1.93 fixed 1.5 x 1032 



"Unabsorbed 0.3-10 keV luminosity at a source distance of 3.5 kpc. 
''Overall reduced x^ = 106 with 1213 degrees of freedom. 
■^Reduced x^ — 0.75 with 16 degrees of freedom. 
''Reduced x^ ~ 1-38 with 7 degrees of freedom. 
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Fig. 3. — Spectral parameters evolution across the outburst. Circled values indicate a fixed param- 
eter in the fitting procedure. 



3. Discussion and conclusions 

SAX J 1808 is the prototype of the accreting millisecond X-ray pulsar (AMP) class. The 
source underwent several outbursts to date which have been monitored in great details thanks to 
RXTE observations (Wijnands 2005; Wijnands et al. 2003). These observations however were 
limited to the brightest part of the outburst, being the PCA on RXTE a coUimated instrument (and 
therefore heavily background limited). Here we report on the first campaign aimed at studying 
the faintest tail of the outburst and return to quiescence. Thanks to the fast repointing and flexible 
scheduling capabilities of the Swift satellite, we monitored the return to quiescence of SAX J1808 
during the 2006 outburst. The bright phase of the outburst is similar to what has been observed in 
the past: following the peak and the smooth decay a flaring behaviour sets in with at least three 
rebrightening episodes. After this the source turns to quiescence even if same low level activity 
(factor of ~ 5) is still present. 

The puzzling result of this observational campaign on SAX J1808 is that out of 13 observa- 
tions with luminosity below ~ 10^^ erg s^^, 8 times the source was found with a luminosity around 
~ 5 X 10^^ erg s^^ (see Fig. 1), i.e. at a luminosity level ~ 10 times higher than the true quiescent 
level (Campana et al. 2002; Heinke et al. 2007). Given the wild luminosity variations this is some- 
what strange and calls for the presence of a 'metastable' state. This luminosity level is difficult to 
interpret in case of standard disk accretion but finds a natural explanation if the accretion process 
is mediated by the magnetic field of the fast spinning neutron star. When the magneto spheric ra- 
dius (i.e. the radius at which the neutron star magnetic field starts controlling the accretion flow) 
is larger than the corotation radius, the system is in the (so-called) propeller regime. In this case 
matter is halted from the rotating magnetosphere and does not accrete onto the neutron star sur- 
face (see e.g. Campana & Stella 2000 and references therein). When the magneto spheric radius 
becomes larger than the light cylinder radius (ric = c P/2 n, where c is the velocity of light and 
P the spin period), the magnetic field cannot corotate any more with the neutron star and a dipole 

3/2 —3/2 

losses will take place. The luminosity range spanned in the propeller regime is ~ 440 P2.5 ^1.4 
(where P2.5 is the spin period in units of 2.5 ms and Mi 4 the neutron star mass in units of 1.4 M©; 
see Campana & Stella 2000). In order to have the 'metastable' state to coincide with the lowest 
luminosity in the propeller regime (i.e. just before the reactivation of the pulsar), we have to re- 
quire a magnetic field of ~ 7 x 10^ G, in line with previous estimates: ~ (1 — 10) x 10^ G using 
disk-magneto sphere interaction models (Psaltis & Chakrabarty 1999), ~ (1 — 6) x 10^ G using 
simple considerations on the position of the magneto spheric radius during quiescent periods (Di 
Salvo & Burden 2003); (0.4 — 1.5) x 10^ G modelling period changes during outbursts with a 
magnetic torque model (Hartman et al. 2008). 

The lowest luminosity level in the proposed interpretation is ascribed to the interaction of a 
turned-on pulsar with the interbinary and circumstellar environment. Variability is expected at this 



10 



stage due to the rapidly changing environment as observed, e.g., in the millisecond radio pulsar 
PSR J1740-5340 (D'Amico et al. 2001; Ferrario et al. 2001) which gets eclipsed over a range of 
orbital phases for different orbital cycles. 
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